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r-> : 

In the framework of QCD sum rules one uses a scheme, allowing 
one to apply the conditions of both nonrelativistic heavy quark motion 
. inside mesons and independence of nonsplitting nS-state density on 

| the heavy quark flavours. In the leading order an analitic expression 

\q ■ is derived for leptonic constants of both heavy quarkonia and heavy 

. mesons with a single heavy quark. The expression allows one explicitly 

^ ' to determine scaling properties of the constants. 

:a; 

Introduction 

By definition, heavy quarks have the ttiq mass values that are much greater 
j> | than the confinement energy A. Therefore, in some cases, consideration of 

hadrons with the heavy quarks allows one to use expansions of some quanti- 
ties over the small parameter of the A/w,q ratio. If the heavy quark virtuali- 
ties inside a hadron are not large, then one allows the kinematical expansion 



vk ~ 0, (1) 
\k 2 \ < m? 



where v is the quark 4-velocity, Pq is the quark momentum, k character- 
izes the heavy quark virtuality. In the kinematics of eq.(|l|) in QCD, the 
heavy quark action, expanded over the small parameter, leads to Effective 
Heavy Quark Theory (EHQT) 0, so, in the leading approximation the the- 
ory possesses the symmetry with respect to the substitution of a heavy quark, 



1 



moving with the velocity v, by any other heavy quark, moving with the same 
velocity v and having an arbitrary orientation of its spin. 

In the system, where v = (1,0), the heavy quark hamiltonian in an ex- 
ternal field has the form 



SB 



where V is the potential, B is the chromomagnetic field. 
In heavy mesons (Qq) with a single heavy quark, one has 



<k 2 > 

m Q 



0(l/m Q ) , (3) 



9 <G B> ~ 0(l/m Q ) , 
m Q 

and the distance r is determined by the light quark motion around the static 
source of the gluon field. Hence, in the leading approximation of EHQT, the 
heavy meson wave function is universal and independent of the flavour of the 
heavy quark inside the meson. This feature leads to both the scaling law for 
the leptonic constants of the heavy mesons 

f . M = const. , (4) 

and the universality of form factors for semileptonic transitions between the 
hadrons, containing a single heavy quark (for example, B — > D^*'lu) [|l|. 

In the case of heavy quarkonium (QQ'), the chromomagnrtic field arises 
only at nonzero velocity of the source, so that 

B ~ 0(v) ~ 0(l/m Q ) , (5) 

and, hence, spin-dependent splittings of the quarkonium levels arise only at 
the second order over I/tuq, so that in what follows, we neglect the spin- 
dependent splittings in the heavy quarkonium. 

As for the kinetic energy of the heavy quark motion inside the quarkonium 

(QQ') 

P P 2 

T = + , (6) 

2m Q 2m Q 
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it determines essentially both the quark binding energy E and the wave func- 
tion ^>E(f), as one can know this fact from an experience of working with 
nonrelativistic potential models of the quarkonia. Moreover, the distance be- 
tween the heavy quarks inside the quarkonium depends on the quark masses. 
Thus, generally speaking, the leading approximation of EHQT may not be 
applied to the heavy quarkonium (QQ'), whose wave function depends on 
the quarkonium content. 

However, as it has been shown in ref.||, in the region of average dis- 
tances between the heavy quarks inside the (cc) charmonium and the (bb) 
bottomonium 

0.1 fm < r < 1 fm , (7) 

and with accuracy up to an additive shift, the QCD-motivated flavour- in- 
dependent heavy quark potentials, behaving as the Coulomb interaction at 
small distances and having linearly rising confining part at large distances 
(Cornell model ||, Richardson potential [§], Buchmiiller-Tye model ||), al- 
low the parameterizations in the forms of logarithmic || and power |7| laws, 
possessing simple scaling properties 

V L {r) = c L + d L \nAr, (8) 
V M (r) = -c M + d M (Ar) k . (9) 

By the virial theorem for average values of the kinetic energies in potentials 
(D, (g), one can get 

< T L > = d L /2 = const. , (10) 
<T M > = j^(c M + E), (11) 

respectively, so that at \E\ <S cm, k <C 1 and with the accuracy by the small 
binding energy of the quarks inside the quarkonium, one concludes 

< T M >^ const. (12) 

In accordance with the Feynman-Hellmann theorem 

(13, 

djj, jj, 
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Table 1: The Mass Difference (in MeV) for the Lightest Vector States with 
The Prescribed Valence Quark Contents. 



state 


T 


i> 


4> 


AM 


563 


588 


660 



where u is the reduced mass of the heavy quarks (QQ f ), and by condition 
([TOD, for the difference of the energies of two levels, one gets 

E{n, fi) - E{n, //) = E{n, //) - E(n, u') , (14) 

i.e. the level density of the (QQ') system does not depend on the heavy 
quark flavours 

dn 

— = const. , (15) 

that is rather accurately confirmed empirically || (see table |ip. 
In the framework of the QCD sum rules || the use of 

1) the small parameter , A/m,Q <C 1, 

2) the nonrelativistic motion of the heavy quarks, v — > 0, 

3) the universality of the quarkonium state density fllSD , 
allows one, in the leading order, 

1) to neglect power A/ibq corrections from quark-gluon condensates, 

2) to take into the account Coulomb-like interactions over ats/v, 

3) to derive the scaling relation for the leptonic decay constants of heavy 

S-wave quarkonium (see ref. [ 



P 

— = const. , (16) 
M K ' 

in the regime, when \niQ — mQi \ is restricted at ?t^q,q' ^> A. Expression 
(|l~6|) is in a good agreement with the experimental values of fr, f\p and 
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In the present paper we generalize the analisys, made in ref.[]TO[, for the 
regime of itlq = xttiqi ^> A and derive the following expression, determining 
the scaling properties of the leptonic constants for the S-wave quarkonia 



2 / h/T \ 2 

= const. , (17) 



M U/J 

where \i is the reduced mass of the quarks. 
Expression (|T7j) 

1) is reduced to eq. (|16D at x — 1, i.e. at 4/i/M ~ 1, 

2) agrees with the scaling law for the leptonic constants of heavy (Qq) 

mesons in the regime M — > oo, /x = const., and 

3) allows one to predict the /b c value for the heavy quarkonium B c with 

the open charm and beauty (the B c search is processed at LEP and 
FNAL). 

In Section 1 the QCD sum rule scheme is considered. It allows one explic- 
itly to use condition (|1^) and the dn/dM n quantity as the phenomenological 
parameter, and to avoid unphysical dependence of results on the external 
parameters such as the number of the spectral density moment or the Borel 
transformation parameter. Expression (|17|) is derived. 

In Section 2 one analyses the scaling relation (|T7j), and in Conclusion the 
obtained results are discussed. 



1 Quarkonium sum rules 

Let us consider the two-point correlator functions of the quark currents 

IV(g 2 ) = tJd i xe t ' lx <0\TJ fl (x)Jl(0)\0>, (18) 
U P {q 2 ) = i J d 4 xe iqx < 0|TJ 5 (x)4(0)|0 > , (19) 



where 



Mx) = Qi{x)^Q 2 {x) , (20) 
Mx) = Qi(x) l5 Q 2 (x) , (21) 

(22) 
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Qi is spinor field of the heavy quark with i = c,b. 
Further, write down 

IV = (~9,u + ^)ny(g 2 ) + ^IW) , (23) 

where Uy and Us are the vector and scalar correlator functions, respectively. 
In the following we will consider the vector and pseudoscalar correlators: 
U v (q 2 ) and U P (q 2 ). 

Define the leptonic constants fy and fp 

<0\J,(x)\V(X)> = ie^fyM v e ik \ (24) 
<0\J 5ll (x)\P> = ik^f P e ikx , (25) 



where 
so that 



Jbp(b) = Qi(xh5l^Q2(x) , (26) 
f M 2 

< 0\J 5 (x)\P >= % JP p e ikx , (27) 
m 1 + m 2 

where \V > and \P > are the state vectors of the 1~ and 0~ quarkonia, 
and A is the vector quarkonium polarization, k is 4-momentum of the meson, 
k 2 — M 2 

npy — 1V1 py. 

Considering the charmonium (ip, ip' ...) and bottomonium (T, T', T" 
...), one can easily show that the relation between the width of the leptonic 
decay V — > e + e~ and fy has the form 

T(V - e + e~) = f e 2 a 2 em & , (28) 

where e$ is the electric charge of the quark i. 

In the region of the narrow nonoverlapping resonances, it follows from 
Eqs.(PI) - (E7D that 



-Smll^V) = J2fvnM v J(q 2 -M Vn ), (29) 



ismlirV) = J2fpnM Pn ( I , 2 % 2 ~M|J, (30) 

n (mi+m 2 ) z 



7T 
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Thus, for the observed spectral function one has 

-3m<f (g 2 ) = -3m< e p s V) + PvMAr) , (31) 

71 71 

where p(q 2 , fi 2 ) is the continuum contribution, which is not equal to zero at 
q 2 > n 2 . 

Moreover, the operator product expansion gives 
Ii^ CD \q 2 ) = n^V) + C G (q 2 ) < —G 2 > +C t (q 2 ) < rmQiQi >+..., 

71 

(32) 

where the perturbative contribution U^ pert \q 2 ) is labled, and the nonpertur- 
bative one is expressed in the form of the quark-gluon condensate sum with 
the Wilson's coefficients, which may be calculated in the QCD perturbative 
theory. 

In Eq.([32]) we were restricted by the contribution of the vacuum expec- 
tation values for the operators with dimension d = 4. For Co'\q 2 ) one has, 
for instance, 

C (P) = 1 <f_ ( 3(3v 2 + l)(l-v 2 ) 2 l+v_ _ 9v 4 + 4v 2 + 3 \ 
G ~~ 192mim 2 q 2 V 2v 5 ^ l-v v 4 J ' 1 ' 

where 

f = g 2 - mi - m 2 2 w 2 = 1 l_i . 34 

T 

The analogous formulae for other Wilson's coefficients can be found in Ref . . 
In the following it will be clear that the explicit form of the coefficients has 
no significant meaning for the present consideration. 

In the leading order of the QCD perturbation theory it was found for the 
imaginary part of the correlator that |9| 

^mlly ert \q 2 ) = - S 2 (3ss — s 2 + Qmim2S — 2m\s)9(s — (mi + m2) 2 |3,5) 



QmnP( 9 2 ) = -^(s~(m 1 -m 2 ) 2 )e(s-(m l +m 2 ) 2 ) , (36) 
8ns z 

where s = s — m\ + to 2 ,, s 2 = s 2 — 4to|s. 

The one- loop contribution into QmU(q 2 ) can be included into the consid- 
eration (see, for example, Ref.p]). However, we note that the more essential 
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correction is that of summing a set over the powers of (a s /v), where v is 
defined in Eq . (j34|) and is a relative quark velocity, and 0:5 is the QCD in- 
teraction constant. In Ref. || it has been shown that the account of the 
coulomb-like gluonic interaction between the quarks leads to the factor 

77/ n 47r 01 s 1 

F(V) = I7 l-exp(-^) ' (3?) 

so that the expansion of the F(v) over as/v <C 1 restores, precisely, the 
one- loop 0(^ & ) correction 

F(«)«l- T ^ ... (38) 

In accordance with the dispersion relation one has the QCD sum rules, which 
state that, in average, it is true that, at least, at q 2 < 



- / ^ds = U {QCD) (q 2 ) , 39 

Ti J s — q 



where the necessary subtractions are omitted. Smn'"^ 2 ) and U^ CD \q 2 ) 
are defined by Eqs.(|2|) - (g|) and Eqs.® - (H), respectively. Eq.(H) is the 
base to develop the sum rule approach in the forms of the correlator function 
moments and of the Borel transform analysis (see Ref.|D[]). The truncation 
of the set in the right hand side of Eq.(p9|) leads to the mentioned unphysical 
dependence of the fpy values on the external parameter of the sum rule 
scheme. 

Further, let us use the conditions, simplifying the consideration due to 
the heavy quarkonium. 

1.1 Nonperturbative Contribution 

We assume that, in the limit of the very heavy quark mass, the power correc- 
tions of the nonperturbative contribution are small. From Eq.([33|) one can 
see that, for example, 

4 P V) « 0{-^—) , A/m lj2 «l, (40) 

mi77l2 
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where v is fixed, q 2 ~ (mi + 1712) 2 , when Qmll^^g 2 ) ~ (mi + mj) 2 . It is 
evident that, due to the purely dimensional consideration, one can believe 
that the Wilson's coefficients tend to zero as \/m\ 2 - 

Thus, the limit of the very large heavy quark mass implies that one can 
neglect the quark-gluon condensate contribution. 



1.2 Nonrelativistic Quark Motion 

The nonrelativistic quark motion implies that, in the resonant region, one 
has, in accordance with Eq.(^), that 

v -»• . (41) 

So, one can easily find that in the leading order 

ZmTir rt \s) « 3mllf rt) ( S ) - {^f , (42) 

so that with account of the coulomb factor 

F(v) ~ ^ ^ , (43) 
3 v 

one obtaines 

^ rt) (^| S g) 2 . (44) 

1.3 "Smooth Average Value" Scheme of the Sum Rules 

As for the hadronic part of the correlator, one can write down for the narrow 
vector resonance contribution 



n^V) = / —T~2 £ fvn M Vn S ( S - M Vn) > (45) 
J s q n 

nirV) = l— 2 Y,fln ( M t ,A S - M Pn): (46) 



The integrals in Eqs.(|45|)-(f46|) are simply calculated, and this procedure is 
generally used. 
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In the presented scheme, let us introduce the function of the state number 
n(s), so that 

n(ml) = k . (47) 

This definition seems to be reasonable in the resonant region. Then one has, 
for example, that 

l - ^mUp s \s) = sf Vn{s) ±J2^-M 2 Vk ) . (48) 
7T w as " 

Further, it is evident that 

^(s-MD^l^ein-U), (49, 
and Eq.(f4"5"|) may be rewritten as 

nrVH/^/^^X;^-*). (50) 

The "smooth average value" scheme means that 

It is evident that, in average, the first derivative of the step-like function in 
the resonant region is equal to 

<^-£*(»- *)>-!• (52) 

k 

Thus, in the scheme one has 
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s — q z 



<nP(? 2 )> - /rf^^w^, (53) 



<n£"V)> - / ^ ^)_ 

J s — q 2 (m 1 + m 2 y as 

Eqs.fl53D-([54]) give the average correlators for the vector and pseudoscalar 
mesons, therefore, due to Eq.(^) we state that 

3m < rt hadr \q 2 ) >= ^mU^ QCD \q 2 ) , (55) 
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that gives with account of Eqs.(|3), (ET) and at the physical points 



M 2 

n 



M n Ti dn \MJ ' ^ 

where in the limit of the heavy quarks we use, that for the lightest resonances 
one has 

m 1 + m 2 ^ M , (57) 

so that 

fvn ^ f'Pn = fn ■ (58) 

Thus, one can conclude that the QCD sum rules give for the heavy quarkonia 
the identity of the fp and fy values for the lightest pseudoscalar and vector 
states. 

Eq.(^) differs from the ordinary sum rule scheme because it does not 
contain the parameters, which are external to QCD. The quantity dM n /dn 
is purely phenomenological. It defines the average mass difference between 
the nearest levels with the identical quantum numbers. 

Note, in the Borel sum rules, the derivative procedure over a gives the 
possibility to find both the constants / and the bound state mass versus the 
current quark mass choice. 

It must be noted, that the approximation made implies that we neglected 
the continuum contribution in the resonant region, and this assumption is 
valid for the lightest states only. 

The relations, connecting / with the dM n /dn value, were derived in some 



other ways, so in the quasiclassical approximation [11], in the sum rule analy- 



sis with the use of the Euler-McLohren transformation [TI| and by the double 



action of the Borel transformation 13 



2 Sum rule analysis 

As it has been noted in Introduction, the phenomenological properties of 
the heavy quark potential lead to the independence of the quarkonium state 
density on the heavy quark flavours (see eq.(|l5|)). Thus, in accordance with 
eq . (|56|) and in the leading approximation with no account of the logarithmic 
and power corrections, one can draw the conclusion, that for the leptonic 
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Table 2: The Experimental Values of the Leptonic Constants (in MeV) for 
the Quarkonia in Comparison with the Estimates of Present Model. 



quantity 


exp. 


present 


f<t> 


232 ±5 


230 ± 25 


u 


409 ± 13 


400 ± 40 


/t 


714 ± 14 


700 ± 70 



constants of the S-wave quarkonia the scaling relation takes place 

p /My 



M \AnJ 



const. , (59) 



independently of the heavy quark flavours, so that the constant in the right 
hand side of eq.([59|) is determined by the expression 

const. = ^ d -Mn. ( 6 0) 
7r an 

In ref. [0], for the numerical estimate of the constant we have used the data 
on the average mass difference of the S-wave bottomonium 

< >= -((M T ' - M r ) + (M T » - M T /)) , (61) 

and the flavour- independent value of as, determining the Coulomb term of 
the Cornel potential, 

a s ~ 0.36 . (62) 

In the case of the quarkonia with the hidden flavour, one has 4/i/M = 1. 
So, the calculated values of the leptonic constants for the T-, ip- and 0- 
mesons are presented in table 0, and they are in a good agreement with the 
experimental values of these quantities. 

For the heavy (be) quarkonium with the open charm and beauty, one can 
easily find the estimate 

f Bc = 460 ± 60 MeV , (63) 
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where the uncertainty is caused by the ambiguity in the choice of the quark 
masses @ 



m c = 1.4 -r 1.8 GeV, (64) 
m b = 4.6^5.2 GeV. (65) 

In the potential models |15, |l^, [T^, |i"8"l , the mass estimates for the S-wave 



levels of the B c mesons agree with the expected behaviuor, when the mass 
difference is practically the same, as for the families of the charmonium and 
the bottomonium. 

The Jb c value (|63|) is in an agreement with the other estimates, obtained 
in the framework of both the potential models |IR IIEL O, II9L EDI, EIL E2 



and the QCD sum rules [[TJ], [16], [18], [23[] , where, in the other schemes, a large 
spread in the Jb c predictions takes place due to the ambiguity in the choice 
of the hadronic continuum threshold, the number of the spectral density 
momentum or the Borel parameter. 

Further, in the limit case of B- and D-mesons, when the heavy quark 
mass is much greater than the light quark mass rriQ ^> m q , one has 

u ~ m q 

and 

f M = — /I . (66) 

7r an 

Then it is evident that at one and the same \x one gets 

f M = const. (67) 



Scaling law (p7|) is very well known in EHQT |IJ for mesons with a single 
heavy quark (Qq), and it follows, for example, from the identity of the B- 
and D-meson wave functions in the limit, when infinitely heavy quark can 
be considered static source of gluon field. 

In our derivation of eqs. (|66D and (^7|) we have neglected power corrections 
over the inverse heavy quark mass. Moreover, we have used the presentation 
about the light constituent quark with the mass, equal to 

m q ~ 330 MeV , (68) 
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so that this quark has to be considered as nonrelativistic one v — > 0, and the 
following conditions take place 

m Q + m q ~ M[qq , m q < m Q , (69) 

and 

fv^fp = f- (70) 
In agreement with eqs. flBTf) and (^), one finds the estimates^ 

f BM = 120 ± 20 MeV , (71) 
/ dW = 220 ± 30 MeV , (72) 

that is in an agreement with the estimates in the other schemes of the QCD 



sum rules 24 



Thus, in the limits of 4/i/M = 1 and /x/M -C 1, scaling law (|59|) is 
consistent. 



In ref. []T4j the sum rule scheme with the double Borel transform has been 
used, so, it allows one to study effects, related to power corrections from 
the gluon condensate, corrections due to nonzero quark velocity and nonzero 
binding energy of the quarks in the quarkonium. 

The numerical effect from the mentioned corrections considers to be not 
large (the power corrections are of the order of 10%), and the uncertainty, 
connected to the choice of the quark masses, dominates in the error of the 
fB c value determination. 



Conclusion 

In the present paper, in the framework of the QCD sum rules and in the lead- 
ing approximation, we have considered the scaling properties of the leptonic 
constants for the S-wave quarkonia with heavy quarks, and in the specific 
scheme, allowing one to use the spectroscopic data on the quarkonium level 
density, we have gotten the relation 

P (M 
M V4^ 

1 In ref. the dependence of the S-wave state density dn/dM n on the reduced mass 
of the system with the Martin potential has been found by the Bohr-Sommerfeld quanti- 
zation, so that at the step from (bb) to (bq), the density changes less than about 15%. 



J = const. , 
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that is in the good agreement with the experimental data on the leptonic 
constants for the T-, ip- and 0- mesons. It allows one to predict the Jb c value 
for the B c meson, whose search is processed at LEP and FNAL. 

The estimates for the leptonic constants of the B and D mesons due 
to the scaling relation is in the agreement with the values, obtained in the 
framework of the other schemes of the QCD sum rules. 
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